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SUMMARY 


The original explanation of shaft whipping due to internal friction by 
A. L. Kimball in 1923 employed a physical argument based on hysteresis in the 
exural stress-strain relation. Ten years later an even simpler physical 
explanation of the destabilizing nature of rotating damping was briefly de- 
scribed by D. M. Smith. Subsequently, there have been many discussions of the 
effects of rotating damping on rotor stability. In most cases, the emphasis 
as been on the mathematical aspects of the stability analysis rather than on 
p ysical explanations. The present note provides an elaboration of Smith's 
physical argument . It gives a clear insight into the destabilizing mechanism 
and permits an elementary physical determination of the stability limit. 


INTRODUCTION 


The fundamental critical speed of an elastic rotor is a kind of natural 
barrier for rotational speed much like the velocity of sound is for an airplane. 
In both cases, there were historical periods when some engineers were of the 
opinion that safe operation would not be possible beyond the barrier. Today 
both barriers are routinely broken but not without careful engineering to 
ensure safety in the supercritical or supersonic regimes. 

The first analytical treatments of supercritical rotor whirling were given 
by Foppl in 1895, for an undamped model, (ref. 1) and by Jeffcoat in 1919, for 
a modei with damping, (ref. 2). In this case, art preceded science, for in 
1895 some commercial centrifuges and steam turbines were already running super- 
critically (ref. 1). These early analyses indicated that, if anything, super- 
critical operation should be smoother than subcritical. Since then, many 
instability mechanisms (ref. 3) have been encountered, each of which poses a 
potential threat to the hoped-for smoothness of supercritical operation. The 
first of these to be diagnosed was the mechanism of internal friction or damping 
of relative motion with respect to the rotating system. The theory was initially 
described in 1923 (ref. 4) by A. L. Kimball, Jr. A fuller account is given in 
reference 5. This is a simple physical argument to show that damping forces in 
the fibers of the bending shaft have a resultant action which builds up a 
whirling motion when the rotation is supercritical. The argument is independent 
of the specific type of damping: viscous, hysteretic, or otherwise. Essentially 
the same argument appears in references 3 and 6. In 1933, Baker included the 
damping action due to linear frequency-independent hysteresis in the equations of 
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motion and determined the growth rate of the unstable whirl (ref. 7). An 
extended presentation of the same analysis appears in reference 8. 

In 1933, there also appeared the remarkable paper of D. M. Smith (ref. 9). 
Many modern investigations in rotor dynamics were anticipated in this paper. 

In particular, instabilities due to unsymmetr ical rotating stiffness and 
unsymmetrical rotating inertia are described for the first time. Instability 
due to rotating damping is analyzed for a symmetrical system and the threshold 
frequency for viscous damping obtained for the first time. An analytical 
explanation for the ameliorating action of anisotropic support stiffness is also 
given for the first time. The basic idea of the physical explanation of de- 
stabilization due to rotating damping which is described in the present note is 
squeezed into a single sentence of Smith 1 s paper at the conclusion of his 
analysis for the threshold frequency marking the transition from stability to 
instability. He also used the idea to give a simple physical explanation for 
the reason why anisotropic support stiffness raises the onset frequency for 
instability due to rotating damping. 

Smith also made a number of shrewd predictions about rotor systems with 
several critical speeds, many of which have been verified. In the case of 
rotating damping, the situation has subsequently been shown (refs. 10,11) to be 
somewhat more complicated than he envisaged. Later investigators have also 
extended the theory of rotating damping to include damping mechanisms with 
arbitrary frequency dependence and nonlinear amplitude dependence (ref. 12). 

At present, almost every engineer dealing with rotor dynamics is aware that 
rotating damping is destabilizing although to many it still is something of a 
mystery that a mechanism of energy dissipation can cause instability. Kimball s 
explanation, while physical and easy to follow step by step, is sufficiently 
intricate because of the awkward geometry of bending during simultaneous rotation 
and whirling that the argument as a whole is not as compelling as it might be. 
Smith’s explanation is much simpler and persuasive but is not widely known in 
America. In 1970, the present author independently developed an explanation 
similar to Smith’s which he presented several times in invited lectures based 
on reference 13. For nearly ten years, he believed his explanation was original 
until on reviewing A. B. Pippard’s recent book (ref. 14) for the Journal of 
Applied Mechanics (ref. 15) he discovered that the explanation was, in fact, 
well known in England and was led back to Smith’s work. 


PLANAR ROTOR MODEL 


Consider the planar model of a rotor mounted on a flexible shaft shown in 
Fig. 1. The rigid ring is forced to turn at the angular rate Q by an external 
source. A mass particle m is suspended elastically from the ring by four mass- 
less linear springs. The mass particle represents the central rotor of the 
classical Foppl-Jef fcoat model and the springs represent the flexible shaft. 
Gravity in the plane of the diagram is neglected. The elastic force system for 
small displacements is circularly isotropic. The equilibrium position is with 
the mass at the origin and when the mass has a radial displacement r, in any 
direction within the plane, the elastic restoring force is kr directed back 


370 




Fig. 1 Planar model of ideal balanced rotor with no damping. 


toward the origin. Furthermore, this characteristic of the elastic system is 
independent of the rotational speed The equations of motion for the mass 

particle in Fig. 1 in terms of the stationary coordinates x,y of Fig. 2 are 


1* I 

1 1 x ) 


m 1 

+ k =0 

(i) 

( y ) 

(y ) 



independently of Q . All natural motions are linear combinations of two inde- 
pendent modes which have the same natural frequency <n n = VkTm. The two 
independent modes may be taken as rectilinear oscillations along two diameters 
or as a pair of circular whirling modes, one counter-clockwise and the other 
clockwise. The resulting orbit, in the general case, is an elliptical path 
which is periodically circumnavigated with the period 27T/ui n . Due to the isotropy 
of the system, it is always possible, without loss of generality, to choose the 
orientation of the x,y axes and the starting phase so that the general free 
motion has the representation 


x = a cos ai n t 


( 2 ) 


y = b sin u t 

where j a | and |b| are the major and minor semiaxes of the ellipse. It is some- 
times convenient to consider the x,y plane as the locus of the complex variable 


z = x + iy ( 3 ) 

In such a setting a counter-c] ockvise whirl of unit amplitude and angular rate to 
is represented by 
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Fig. 2 Stationary coordinates x,y and rotating coordinates s,T] f° r 
the mass particle in Fig. 1. 


e 1(AJt = cos cot + i sin cot (^) 

The elliptical motion (2) can be represented as a superposition of counter 
clockwise and clockwise whirls by inserting ( 2 ) into (3) and expressing the 
trigonometric functions in terms of exponentials of imaginary argument. The 
result is 


z 


) s co^t + ib sin oo n t 


. icOr-vt , -pj iw n t 

A e 11 + Be 11 


(5) 


where 


A = 1/2 (a + b) a = A + B 

B = l/2(a - b) b = A - B 


( 6 ) 


While not necessary for the physical explanation which follows, the equations of 
free motion with respect to the rotating axes £,n of Fig. 2 are given here for 
completeness 


m 





hi 


*+■ 2m^ 


0 

1 



(7) 


In the elliptical orbit (2), the instantaneous kinetic energy is 

T = S’"" " 4 - b " cos 2u) nt) W 

and the intantaneous potential energy is 

v = k ( a 4 — + 4 — cos 2 w n c ) (9) 

These energies oscillate about the same mean value with opposite phase so that 
the total energy E in the rotor orbit remains constant 
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E = T + V = l/2k(a 2 + b 2 ) = k (A 2 + B 2 ) 


( 10 ) 


QUALITATIVE EFFECTS OF DAMPING 


Consider now the introduction of isotropic linear damping to the rotor 
model of Fig. 1. In the following section, arbitrary frequency dependent 
linear damping mechanisms will be treated. For the present, it is sufficient to 
consider classical viscous damping. The damping of the motion of m may be with 
respect to the stationary axes x,y or the rotating axes £,n in Fig. 2. Damping 
with respect to the stationary axes is represented in Fig. 3 by the four dashpots 
whose outer extremities are fixed in the stationary reference frame. For small 
motions of m, the system of dashpots develops a circularly isotropic damping 
force. If the absolute velocity of m is v in any direction within the plane of 
motion, the resultant dashpot force is c s v in the opposite direction. In 
particular, if m is forced to travel in a circular orbit of radius A with counter- 
clockwise angular rate co n , the damping force c s AoJ n acts tangent to the circle in 
a clockwise sense. 

Damping of relative motion with respect to the rotating system is represented 
in Fig. 4 by the four dashpots whose outer extremities are fixed in the rotating 
ring. For small motions of m, this set of dashpots develops a circularly 
isotropic damping force which acts in opposition to the velocity of m relative 
to the rotating system. In particular, if m is forced to travel in a circular 
orbit of radius A with counter-clockwise angular rate co n as viewed from the 



Fig. 3 Damping of motion with respect to stationary axes is 
represented by four symmetrically oriented dashpots. 
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stationary axes, it will appear to have a counter-clockwise angular rate tn n U 
in the rotating frame. The damping force generated by the rotating dashpots is 
c r A(w n - ft) acting tangentially to the circle in a clockwise sense. 

The effects of these damping mechanisms on the free motion of the rotor can 
be argued qualitatively as follows. Suppose that the damping forces are very 
small in comparison with the spring forces. This implies that a free orbit will 
undergo only a small change during one period 2ir/u) n as a result of the damping. 

If initially the orbit is a counter-clockwise circle of radius A, the preceding 
discussion shows that when the rotor rotation rate SI is subcritical; i.e., less 
than the natural frequency of the free motion (n n , both the stationary and rotating 
dashpots act to retard the motion of m. The rotor does work on the dashpots and 
the total energy E of the rotor orbit is diminished. If initially the orbit is 
a superposition of a counter-clockwise whirl of radius A and a clockwise whirl 
of radius B as in (5) a similar argument shows that when ft < w n the rotor does 
work on the dashpots in both whirls and that both A and B are diminished after 
a period 2TT/u) n . Any rotor orbit generated by an accidental disturbance will thus 
be damped out and the system is stable. 

When the system rotates supercritically , the rotating dashpots do work on 
the rotor and add energy to the rotor orbit. This destabilizing action can be 
seen by returning to Fig. 4 and considering the case of a counter-clockwise whirl 
of radius A at absolute angular rate u> n when SI > o) n * The relative motion now is 



Fig. 4 Damping of relative motion with respect to rotating 
axes is represented by four symmetrically oriented 
dashpots . 

a backward (i.e., clockwise) whirl with rate fi - w n . The resultant dashpot 
force acts tangentially to the circular path in a forward (i.e. , counter- 
clockwise) sense with magnitude c r A(ft - io n ) . The supercritical rotation of the 
system thus acts to drag the rotor forward around its orbit through the rotating 
dashpots . 
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j , _ . aC ^° n "ay be seen m an even clearer fashion if, instead of the four 

as po s m Fig. 4, it is imagined that the interior of the ring is filled with 
a massless viscous fluid which rotates with the ring and acts to retard the 
relative motion of the mass m with respect to the ring. When the mass has a 
circular orbit m the same absolute sense as the ring rotation, but at a slower 
rate, the viscous drag pulls the mass forward and adds energy to its orbit. 

The model described by Pippard (ref. 14) is essentially equivalent. The 
r 0t ?j aS modelled ^ a conical pendulum with a heavy bob suspended in the gravity 
leld by a string of length L. The natural modes are forward and backward 

Damnlni ? ! naCural fre ^cy w n = A/g/L (for small cone angles). 

P u 8 modelled hy allowing the pendulum bob to dip into a glass of water. 

the glass is on a variable speed turntable, damping with respect to a rotating 
frame can be demonstrated. "If the fluid is caused to rotate more slowly than 
the pendulum bob in its circular orbit, the bob experiences a retarding force 
and sinks toward the vertical. But if the fluid rotates faster than the bob 

it urges the bob onwards and causes its orbit radius to increase" (ref 14) 

See Fig. 6 below. v * J * 

Returning to the planar rotor model with both stationary damping Fig 3 

u" d '° ta J 1 J g K da, ”P ln8, Fig - 4 ’ a forward whirl at absolute angular rate oo n will 
e retarded by the stationary damping force and will be urged on by the rotating 
damping force when the rotation is supercritical, Q > o) n . Whether the net effect 
is to remove or to add energy to the whirl orbit, depends on which force is 
arger.^ Neutral stability occurs when the forces are equal. There is then no 
change in the orbit energy and the free motion persits indefinitely The 
quantitative determination of the stability borderline is discussed in the 
following section, taking account of the frequency dependence of the damping 


In the case of a backward whirl, both the rotating dashpots and the 
stationary dash P ots act to retard the whirl. The orbit energy of a backward 
whirl is always decreased by either rotating or stationary damping. For 
completeness the equations of motion for the planar rotor model of Fig. 2 with 
the damping elements of Figs. 3 and 4 are given here. In the stationary 
coordinates, the equations for the damped system are 


m 


( x ) ( x ) ( x ) 


0 

-1 


H 

M 


(ii) 


in place of (1) for the undamped system. In the rotating coordinates used in 
w J , the equations for the same damped system are 




+ <c s + c r ) j. j + (k - "~ 2 > | 


+ 2mft 


-1 

0 


In) 


+ c Q 
s 


-1 

0 




'n 


'= 0 ( 12 ) 


Note that in the stationary coordinates of (11), the rotating damping is 
represented by both a normal damping term and by a quasi-gyroscopic term. In 
fashion 3 coordinates (12), the stationary damping is represented in like 
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QUANTITATIVE EFFECTS OF DAMPING 


In the preceding section, it was shown that at the stability borderline for 
a forward whirl the resultant forward drag force of the supercritically rotating 
damping balances the backward drag of the stationary damping. To give a 
quantitative expression of this we introduce the concept of loss facto r. When 
a material or mechanical element undergoing harmonic deformation at frequency 
oj has both elastic and damping behavior, the ratio of the amplitude of the 
damping force to the amplitude of the elastic force is called the loss tangent, 
or loss factor, n. In general, U depends on both the amplitude and frequency 
of the deformation. For linear damping mechanisms, n is independent of amplitude. 
There are many real damping mechanisms that can be satisfactorily modeled as 
linear frequency-dependent processes; e.g., polymers, and metals when the 
dominant damping is due to transverse thermal currents (refs. 13,16). For the 
planar rotor model of Fig. 1, we can define a stationary loss factor U s (id) by 
taking the ratio of the tangential drag force, due to stationary damping, to the 
radial spring force when the mass particle executes a circular orbit with radius 
A at frequency w. When OJ = w n , the backward drag force D s due to the stationary 

damping is 

D = rj (03 ) kA (13) 

s s n 

Similarly, a rotating loss factor n r ( w ) is defined as the ratio of the tangential 
drag force, due to rotating damping, to the radial spring force when the mass 
particle executes a circular orbit with radius A at a circular frequency to, wit _ 
respect to the rotating sy stem . In a forward whirl at absolue rate to n > the 
circular frequency with respect to the system rotating at rate 8 is to n - Thus, 

when the system rotates supercritically, the drag force D r due to the rotating 
damping is forward and has the magnitude 

D = n - ia> ) kA (1^) 

r r n 

for an orbit of radius A. If D r = D s , we have a steady orbit of fixed radius A. 

If D <D, there is net retardation and the orbit energy decreases. If 
D >V there is net acceleration and the orbit energy increases. Alternatively, 
the stability is decided by the relative magnitudes of the rotating and stationary 
loss factors as indicated in Fig. 5. The stability borderline condition occurs 
when the loss factors are equal 

r| (8 - (jJ ) = U (o ) (15) 

r n s n 

The running speeds 8 which satisfy (15) depend crucially on the frequency 
dependence of the rotating damping mechanism and on the magnitude of the stationar 

damping . 

The planar rotor model also provides a clear insight into the energy relation 
ships during whirling of a rotor with stationary and rotating damping mechanisms. 
Consider a forward whirl with radius A and angular rate ta n . If the system is 
rotating supercritially at rate 8, the external agent responsible for maintaining 
the rotation rate 8 delivers energy to the system. This power input is mostly 
dissipated in the stationary and rotating damping mechanisms. If there is excess 
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input power, energy will be added to the whirl orbit, if the input power is 

inadequate, energy will be removed from the whirl orbit. Quantitatively, we have 

the following powers, when Q > oo n 

P. = D Afi = kA 2 ftn 

input r r 


P , . „ = D A(fi - 0 ) ) = kA 2 (n - u )n 

diss,rot r n nr 


P , . = D Ago 

diss,sta s n 


(16) 


= kA go n 
n s 


P 


orbit 


(D-D ) Ao 3 
r s n 


kA 2 w (r| - r\ ) 

nr s 



Fig. 5 Forward whirl is unstable for supercritical rotation speeds £2 at 

which the loss factor of rotating damping is greater than the loss 
factor of stationary damping. 


Note that these satisfy the requirement of energy conservation 

P , =P +P +p (17) 

input diss , rot diss,sta orbit v ' 

At a stability borderline, we have no power flow into, or out of, the orbit which 
leads back to the requirement (15). Note also that if the frequency dependence 
of n r ( w ) is such that the instability onset rotation is very close to 0) n then 
the relative amount of energy dissipated in the rotating damping mechanism will 
be very small in comparison to the energy dissipated in the stationary damping 
mechanism even though it is the rotating damping force which is responsible for 
the instability. 

The power flow into the orbit in (16) can be used to estimate the growth 
rate of the orbit in the unstable regime. Assuming exponential growth of orbit 
radius 
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A = A e 
o 


(18) 


at 

with small growth rate a, the energy added to an orbit with initial radius A q 
and angular rate w n would be approximately 


3 27TkA 2 (n 

orbit w n or 


n s ) 


(19) 


During this period, the radius increases from A q to A q exp {2TTa/o) n } and the 
increase in orbit energy as given by (10) is 


4iTq 

AE = kA 2 (e - 1) = 4ukA 2 — (20) 

o v o w n 

Equating (19) and (20), we find 


a = i/2w n (n r - n s ) (21) 

as the estimated growth rate. In (21), n r (w) is to be evaluated at to = ° - W s 
and r| s ((jo) is to be evaluated at OJ = oo n . 

The proceding results can, of course, also be derived by solving (11) or 
(12). If in (11), for example, we set z = x + iy and = k/m the complex 
variable z(t) must satisfy the equation 

x + £ + (a) 2 _ i£r Q)z = 0 (22) 

m n m 

Exponential terms of the form 

z = A exp { (a + i 3) t } (23) 


satisfy (22) if 


8 2 = < - « 2 


a 


(SI - 3)c r - 8c q 

2P> “ 


(24) 


These relations apply with constant c s and c r if the damping is viscous. For 
other linear frequency-dependent damping mechanisms with loss factors Hs( ! ) an< l 
’ 1 r (u)) we can use the relations (see ref. 13) 

c = ^s(B) _ jEjDx (R ~_.B)_ (25) 

c s | 3| r " 6| K } 

on the assumption that a < < 6. This assumption is verified a posteriori for 
light damping. Then, correct to first order in the loss factors, we have from 
(24) 
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a i 88 l/2w n [ri r (fi “ w n ) sgn (fl - 0 J n ) - n s (^n)] a 2 = -l/2a) n [n r (S + a> n ) + n s (u) n )] (26) 


To this order of accuracy, the general solution to (11) is 


z(t) = A l exp{(a x + iaj n t)} 4* A 2 exp { (a 2 - ito n t)} (27) 

This is a superposition of Archimedean spirals. The forward spiral has a growth 
rate a 1 which is positive if fi > a) n and r| r (fi - co n ) > r} s (u n ) . The magnitude of 
Oil in this case agrees with the estimate (21). If ft < 0 J n then a x is negative 
and the forward spiral decays. The growth rate ot 2 of the backward spiral is 
always negative. 


CONCLUSIONS 


Visualization of damping as a drag force on the orbit provides a clear 
explanation of the destabilizing effect of rotating damping. When the rotation 
is faster than the whirl, rotating damping drags the orbiting particle forward. 
When stationary damping is also present, the stability borderline is readily 
determined by balancing the backward and forward drags. A key notion here is 
that a forward whirl at rate oo n with respect to stationary axes appears to be a 
backward whirl at rate 0, - co n with respect to a system rotating supercr itically 
at rate The growth rate of unstable whirls (or the decay rate of stable 

whirls) is readily estimated by a simple energy balance. 

As a final note, we call attention to D. M. Smith’s extension of the previous 
argument to rotor systems with anisotropic non-rotating support elasticity. 

If the bearings in Fig. 1 are not rigidly supported, but are elastically 
supported so that the compliance in the x - direction is greater (say) than that 
in the y - direction, the free motion of the undamped system becomes a super- 
position of a horizontal oscillation at frequency oo x and a vertical oscillation 
at frequency (j 0 2 > tOj . The more anisotropic the supports, the greater the 
separation between and 0) 2 . The resultant free orbits no longer whirl with a 

fixed sense and "there is no tendency to set up a whirl of the type which 

can be dragged forward by rotary damping until the rotary damping forces have 
been so far increased by rising speed that they are commensurate with the 
difference between elastic restoring forces in the two principal directions 

" (ref. 9). Thus, anisotropic bearing supports should raise the onset speed 

Q for instability due to rotating damping. Smith (ref. 9) and others (ref. 17) 
have verified this by applying Routh’s criterion to the equation of motion. 
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APPENDIX 


The description of the destabilizing action of rotating drag forces given 
above can be extended to give a heuristic explanation of the destabilizing 
tendency of fluid filled journal bearings when the rotational speed exceeds 
twice the natural frequency of the rotor. Here, instead of modelling the rotor 
as a particle travelling in an orbit within a fluid bath whose drag either 
retards or advances the particle, we model the rotor as a cylindrical journal 
of radius r rotating within a cylindrical bearing with a small radial clearance 
h(h « r). When the journal is centered and rotating at angular rate ft, the 
incompressible viscous fluid in the gap is pumped circumferentially (longitudinal 
flow is neglected). The fluid velocity varies linearly across the channel of 
thickness h, from v = rft at the journal to v = 0 at the bearing. This flow 
can be decomposed into a mean flow with uniform velocity v = rft/2 and no 
vorticity and a residual flow with no mean flow and large vorticity. 

Next, consider a small circular whirl of the rotor with amplitude A (A << h) 
and angular rate co n where as before 0J n is the rotor natural frequency. As the 
center of the journal traverses a circle of radius A, a wave of thickness 
variation moves around the channel between the journal and the bearing. At the 
thinnest point, the thickness is h - A and at the thickest point the thickness is 
h + A. If we follow the position of, say, the thinnest point, we note that it 
appears to whirl at the rate u> n ; i.e., the circumferential velocity of the wave 
of thickness variation is ru) n . 

Finally, we conclude that such a natural whirling motion will be retarded 
or aided by the mean flow of the pumped fluid depending on whether the mean 
fluid velocity rft../ 2 is less or greater than the circumferential velocity of the 
wave of thickness variation roJ n . Neglecting other mechanisms, the stability 
borderline is determined by a balance between the destablizing tendency of the 
mean flow and the stabilizing tendencies of the residual high-vor ticity flow. 

A necessary condition for the onset of instability is ft > 2u) n . 

The action of the mean fluid flow in a fluid-filled bearing, in dragging 
around a wave of thickness variation when ft > 2u) n is thus proposed as a 
heuristic physical explanation for the classical oil— whip phenomenon. 
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